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ABSTRACT 

Loops of magnetic field in the corona are observed to oscillate and these oscil- 
lations have been posited to be the superposition of resonant kink waves, trapped 
between the loop's footpoints in the photosphere. To date, most analyses of these 
oscillations have concentrated on calculating the frequency shifts that result from 
spatial variation in the kink wave speed (produced primarily by stratification in 
the density and magnetic field strength). Further, most have ignored gravity 
and treated the loop as a straight tube. Here we ignore spatial variation in the 
wave speed, but self-consistently include the effects of gravity and loop curvature 
in both the equilibrium loop model and in the wave equation that governs the 
propagation of the kink waves that live upon the loop. We model a coronal loop 
as an isolated, thin, magnetic fibril that is anchored at two points in the photo- 
sphere. The equilibrium shape of the loop is determined by a balance between 
magnetic buoyancy and magnetic tension, which is characterized by a Magnetic 
Bond Number e that is typically small |e| << 1. This balance produces a loop 
that has a variable radius of curvature, with the legs being relatively straight 
and the apex of the loop the most curved. Further, a loop becomes unstable to 
buoyant rise if the footpoint separation becomes larger than 2tt times the coro- 
nal pressure scale height. The resonant kink waves of such a loop come in two 
polarizations that are decoupled from each other: waves with motion completely 
within the plane of the loop (normal oscillations or "vertical modes" ) and waves 
with motions that are completely horizontal, perpendicular to the plane of the 
loop (binormal oscillations or "horizonal modes"). We solve for the eigensolu- 
tions of both polarizations using perturbation theory for small Magnetic Bond 
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Number. For modes of the same order, normal oscillations have smaller eigen- 
frequencies than binormal oscillations. The additional forces of buoyancy and 
magnetic tension from the curvature of the loop increase and decrease the mode 
frequencies, respectively. The ratio of the frequencies of the first overtone to 
the fundamental mode — a common diagnostic used to assess the stratification 
of the wave speed along the loop — is modified by the inclusion of buoyancy and 
curvature. We find that the normal polarization possesses a frequency ratio that 
exceeds the canonical value of 2, whereas the binormal polarization has a ratio 
less than 2. 

Subject headings: MHD — waves — Sun: Corona — Sun: magnetic fields - 
Sun: oscillations 



1. Introduction 



I nitially through T RACE and EIT observations, and now ongoing with AIA imagery 



(.e.g. 



White et al.l 120121 ). the loops of magnetic field that often overli e magnetic active re 



gions are seen to oscil l ate with frequencies of typica lly 2-4 mHz (e.g.. ISchrijver et al.l Il999 



Verwichte et al.l l2004t Ivan Doorsselaere et al.l 120071 ). The attribution of these oscillations 
to the modes of a one-dimensional cavity, and the observational identification of overtone 
frequencies of this cavity, has given birth to the field of coronal seismology (see Nakariakov 
& Verwichte 2005 and the references therein). The promise offered by observable resonant 
oscillation frequencies has raised hopes for determining physical parameters of coronal loops 
which cannnot otherwise be measured directly. 

Although multiple claims have been made for t he observation of sausage way es in 



post-flare coronal loops (e.g., ISrivastava et al.l 120081 : iNakariakov fc Verwichtel 120051 ) . the 



majority of the observed oscillations useful for seismological purposes are kink o scillations 
rtSchriiver et~allll999l . I2OO2I lAschwanden et al.l Il9~9~9l . I2OO2I : IVerwichte et allbooJ . The dis- 
persion in the observed overtones, compared to theor etical predictions, have been used mainly 



to estimate the mean magnetic field strength (e.g.. INakariakov fc Ofmanll200ll ). it's longi- 



tudinal gradient ( jVerth fc Erdelvill2008l: iRuderman et al 



inside and outside the loop (iDonnellv et al. 



Orzaet al. 2012; van Doorsselaere et al. 



2006 



2007 



2008 



Diaz et al- 



ii) , or the density stratification 



2004 



Goossens et al 



2007: McEwan et al. 2008 



2006). Observational efforts 



are now focused on improving the measurement of oscillation eigenfrequencies by reducing 
the observational errors and by analy zing additiona l loops . A relevant overview of observa- 
tional measurements can be found in lAndries et al.l (120091 ). 
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On the theoretical side, advancements are steadily being made by including additional 
physical effects with the aim of developing a general model of coro nal loops and their oscilla- 
tions. Most of this work is based on the pioneering calculations of lEdwin fc Roberts! ( 119831 ). 
which were presciently performed prior to the first observations of coronal loop waves. This 
wor k revealed the MH D w ave modes of a trans lationally invariant magnetic cylinder. Since 
then lDiaz et al.l (120041 ) and lAndries et al.l (120051 ) have investigated modifications to the eigen- 
frequencies that arise from density stratification along the loop. The expected dispersion has 
been compared to measured values with the goal of constrainin g coronal loop models (see also 
Goossens et al.ll2006l ; lAndries et al.ll2005tlMcEwan et al.l 20081). Fur ther, the quantitative ef - 
fects of loop geometry were investigated by Dymova et al. ( 2006 ) and Verth &: Erdelvi (I2OO8 ) 



with p articular focus on the curvature of the loop axis. For example, Ivan Doorsselaere et al. 



( 120041 ) have shown that the frequency shift induced by curvature appears at second order 
in the ratio of tube radius to the loop length, which is believed to be quite small for the 
majority of loops. They and others have also sugg ested that since the eigenfrequencies for 
both horizontal and vertical ( jWang fc Solarl 120041 ) polarizations are identical, there is n o 
distinction between the oscillations with different po larization (e.g.. iTerradas et al.l 120061 ). 
unless of course the tube is not locally axisymmetric (IRudermanl 120031 ) . 



Here, our goal is to self-consistently include the effects of gravity and curvature, both 
in the establishment of the static equilibrium of a coronal loop and in the propagation of 
kink waves along the loop. We accomplish both of these goals by exploiting the thin-tube 
approximation, whereby lateral variation across the tube is ignored. Under this assump- 
tion the equilibrium shape of the coronal loop is established by a balance between the 
forces of magnetic buoyancy and magnetic tension. The eigenmode problem is reduced to a 
one-dimensional wave equation which we solve through perturbation theory. The paper is 
organized as follows: in Section 2 we present the equation of motion appropriate for a thin 
magnetic fibril. In Section 3 we derive the equilibrium height and shape of the loop, while 
Section 4 details the solution to the kink wave eigenmode problem. In Section 5, we discuss 
the significance of our findings and the implications of our base assumptions. Finally, we 
state our findings and conclusions in Section 6. 



2. Equation of Motion 

We will treat a coronal loop as a thin magnetic fibril embedded within a field-free coronal 
atmosphere. We recognize that both of these assumptions are suspect. Even though the thin 
tube approximation may only apply to a subset of coronal loops, it provides mathematical 
tractability. The assumption that the loop is an isolated magnetic structure is probably 
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never strictly valid, as the entirety of the solar corona is magnetically dominated and most 
loops are an integral part of larger magnetic structures. In fact, most loops are probably 
just selectively illuminated field lines amongst many that form magnetic arcades. Despite 
these objections, we adopt an isolated fibril model for convenience. The consequences of 
these assumptions are discussed in more detail in Section [5j 

The forc es acting on an isolated, thin magnetic fibril have been previously derived by 



Spruitl (119811 ) through averaging the MHD momentum equation over the cross-sectional area 
of the tube. In terms of the acceleration of the tube, the resulting equation of motion is 
given by 



s + c 2 k + D (s x g) x s , (2.1) 

where all fluid properties represent cross-sectional averages and s and k are the instantaneous 
tangent unit vector and curvature vector of the tube's axis. The gravitational acceleration 
is g and p is the gas pressure. The longitudinal forces (those in the s direction) are just 
the pressure and gravitational forces that act parallel to the magnetic field. The last two 
terms on the right-hand side are the transverse forces of magnetic tension and buoyancy 
with c being the kink wave speed and D the tube's fractional overdensity compared to its 
surroundings, 



dv 
~dt 



P 



-Vp + g 



D = 



B 2 
An (p + p e 

P- Pe 
P+ Pe' 



(2.2) 

(2.3) 



The densities p and p e are, respectively, the mass density internal and external to the tube 
and B is the axial magnetic field strength of the fibril. 

The transverse forces (tension and buoyancy) in equation flZTJ include the effects of 
"enhanced inertia," which represents the backreaction caused by the external medium when 
transverse motions of the tube move surroundi ng fluid out o f the way. The form of the 
enhanced inertia is the same as propounded by ISpruitl (1198 lh . whereby the density in the 
inertial term is simply augmen t ed, p — > p+p P . T his form has inspired much controversy (e.g., 
Choudhuril Il990t IChengi Il992l ; iFan et al.l I1994J ). but faithfully reproduces the proper kink 
wave speed and is consistent with all alternate formulations as long as (1) the background 
atmosphere lacks rotational shear and magnetic field, (2) flow along the tube is not permitted, 
and (3) correction terms that are quadratic in the fluid velocity are ignored. 



-5 - 



3. The Equilibrium Shape of the Loop 



We neglect the spherical geometry of the solar atmosphere and assume that the corona 
can be treated as plane-parallel with constant gravity g. We employ a Cartesian coordinate 
system, with the x-y plane corresponding to the photosphere and the z coordinate increasing 
upwards (i.e., g = —gz). The equilibrium position of the coronal loop is confined to the 
x—z plane with the endpoints of the loop anchored in the photosphere at the coordinates 
(±X,0,0). The height of the loop above the photosphere is given by the function Zq(x), such 
that the loop's axis is traced by the position vector Vq(x) = xx + zo{x)z. Figure [1] provides 
an illustration of the geometry of the loop and the coordinate system. 

In addition to this Cartesian coordinate system, we define a local set of Frenet coor- 
dinates for the equilibrium position of the loop. Let s denote the arclength along the loop 
measured from the footpoint located at x — —X. The longitudinal unit vector, or the tan- 
gent vector that lies parallel to the loop's axis pointing in the direction of increasing s, is 
denoted s . Within the x—z plane and everywhere perpendicular to s is the principle nor- 
mal k . We shall soon see that for the loop model we adopt, only convex loops are stable. 
Therefore, the curvature vector ko always points inward and downward in the direction of 
the principle normal ko, with a magnitude equal to the reciprocal of the local radius of 
curvature Rq 1 of the loop's axis. The third orthogonal direction is given by the binormal 
unit vector b = s x k , which is everywhere constant b = y. In what follows we will need 
the standard geometrical equations which describe the Frenet unit vectors, Sq and ko, and 
the radius of curvature Ro in terms of the loop height, 



dr, 



o 




(3.4) 



ds 



k 



ds, 



ds 








(3.5) 




(3.6) 




(3.7) 



where primes denote differentiation with respect to the photospheric coordinate x. 



Since a loop in a state of static equilibrium is confined to the x-z plane, none of the 
forces on the fibril have a component in the binormal y-direction and we only need to consider 
two components of the equation of motion (12. ip . 
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1 dp . . 
— a ^- s o = 

Po 

c 2 

JL-gD z-k = 



(3.8) 
(3.9) 



In these two equations the subscript indicates a background, equilibrium quantity within 
the loop. The triple cross-product appearing in the buoyancy term in equation ( 12. X[) has been 
reduced to a simpler form by making use of the identities (s x z) x s = z — (z ■ s )s = 
(z ■ k )k . 

Equation ( 13. 8 p expresses the balance of forces in the tangential or axial direction s ? 
whereas equation ( 13.91) describes the balance in the transverse direction of the principle 
normal fco- The axial equation is simply hydrostatic balance along magnetic field lines 
which we will satisfy trivially by assuming that the loop and the surrounding corona are 
both isothermal. Furthermore, since the tube is thin and thermal diffusion wipes out lateral 
variations, the temperatures inside and outside the tube are identical, T —T e . All densities 
and pressures (both gas and magnetic) therefore vary exponentially with height z with a 
common scale height H = R gas To/g- A consequence of this thermal uniformity is that the 
kink speed Cq, the plasma-parameter = 8ttpo/Bq, and the overdensity Dq are all constants. 

The shape of the loop Zq(x) is constrained by the balance of the transverse forces of 
magnetic tension and magnetic buoyancy, which is quantified by equation (13. 9p . Since the 
temperatures inside and outside the loop are identical, and the external fluid is field free, 
the tube must be partially evacuated (Do < 0) in order to ensure pressure continuity across 
the lateral boundary of the flux tube. Thus, magnetic buoyancy tries to lift the loop higher 
into the atmosphere and magnetic tension arising from the convex curvature of the loop tries 
to hold it down. Since the kink speed Cq is constant for our isothermal model, all variation 
along the loop in the magnitude of the magnetic tension force comes from changes in the 
radius of curvature Rq- Similarly, the buoyancy only varies because the transverse direction 
k rotates relative to the vertical z. 

If we insert equations (13. 5p and (13. 6p into equation (I3.9p . we obtain a nonlinear ODE 
for the height of the loop zo(x), 



z'' e 



1 + (z'oY X 
where we have defined a Magnetic Bond Number, 



(3.10) 
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e=— ^. (3.11) 

c o 

The Magnetic Bond Number can be written in a more recognizable form if we multiply the 
numerator and denominator by the enhanced inertia p + p e , 

_ g(p- Pe )X 

In this expression the surface tension that normally appears in the denominator of the Bond 
Number has been replaced by magnetic tension. The Magnetic Bond Number e is a signed, 
dimensionless number that embodies the relative importance of the buoyancy and magnetic 
forces. The Magnetic Bond Number appears in calculations of the Rayleigh- Taylor instability 
when one of the fluid layers is filled with a h orizontal field, providin g the critical wavenumber 



below which instability ensues, k c = —e/X (I Chandr asekharl 1 1 96 ll ) . 



For our isothermal model, e is a constant and furthermore it can be expressed as a 
ratio of length scales, in particular the ratio of the footpoint separation to the pressure 
scale height e = —X/2H. This property is easily derived by realizing that the combination 
of the perfect gas law and pressure continuity across the lateral surface of the tube allows 
one to express the magnetic pressure inside the tube as a ratio of the density difference to 
the fluid temperature. When this expression is inserted into equation (13.121) . the density 
difference cancels in the numerator and denominator, leaving just the pressure scale height 
in the denominator (which is proportional to the temperature). 

Since e is a constant, equation (13.101) can be integrated directly to obtain the derivative 
of the height function z' (x). Similarly, this expression for the derivative can be subsequently 
integrated to obtain the height function itself z (x). We choose the constants of integration 
such that z (±X) = and z' Q (0) = 0. A summary of the resulting properties of the loop is 
provided below: 



z (x) = — ln( C ,° S e , , J , (3.13) 

e \cos (ex/A J / 

z' (x) = tan (ex /X) , (3.14) 



z»(x) = — sec 2 (ex/X) , (3.15) 
X 
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Ro(x) = -—sec (ex/A) , (3.16) 
s(x) - X \n( 1 + Sin£ 1±™&PQ \ (317) 

The last of these equations was obtained by integrating equation (13. 7p . with the constant 
of integration chosen such that s(—X) = 0. Figure [2] displays solutions for several different 
values of e. For reference, the dashed curve corresponds to a semi-circle, clearly illustrating 
the variable radius of curvature for the loops. 

The length of the loop L is obtained by inserting the footpoint position x = X into 
equation ( I3.17p . 

, X ^ /1 + smA _ 



e \ 1 — sin e _ 

Two interesting limits of this equation exist. As e — > the length of the loop converges to 
the footpoint separation L — > IX. This arises because the loop becomes straight, flat, and 
confined to the photosphere. As e — )■ ±7r/2 the loop length diverges logarithmically because 
the height of the loop grows without bound. This can be seen by recognizing that the loop 
reaches its apex at its center [x = 0), therefore achieving a maximum height 



^apex = ^o(O) = — In (cose) . (3.19) 

One can easily see that e = ±7r/2 corresponds to a logarithmic singularity in the height. 

Clearly the height Zq(x) should be a positive function for the range x G (—A, A). With 
a little thought, from equation ( 13. 13[) we can see that this is only possible if e G (— 7r/2, 0). 
Therefore, for constant Co only underdense tubes (D < 0) form stable loops and those 
loops are comprised of a single arch without inflection points (z'q < 0). Furthermore, for zq 
to remain positive, e must be bounded and an equilibrium is only possible if the magnetic 
tension exceeds a threshhold, 



%>% = 2K- 1 g( Pe -p)X. (3.20) 

47T Att 

This condition for stability can also be used to constrain the footpoint separation. Since, 
e = —X/2H, stability requires that the footpoint separation be less than a critical value that 
depends on the corona's pressure scale height, 
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2X < 2X C = 2ttH . (3.21) 

Typically, we do not need to be overly concerned about the limit of large buoyancy forces. 
For most coronal loops the Magnetic Bond Number e is quite small. Using t ypical values of 



the f ootpoint separation 2X = 150 Mm, kink wave speed cq = 1 Mm s 1 (ITomczyk et al. 



20071 ). overdensity D$ = — 1 (appropriate for very small plasma parameter /3), and gravity 
g = 2.7 x 10~ 4 Mm s -2 , we find from equation (13. lip that e = —2 x 10 -2 . For small values 
of e, the loop is short and low-lying with z apex = |e| X/2. This geometry is realized because 
the field lines exit the photosphere at an oblique angle — i.e., z' (±X) = tane. As e becomes 
small, the angle between the fibril's axis and the photospheric surface approaches zero (see 
Figure |2]). 



4. Kink Wave Oscillations 

Kink waves cause the axis of the fibril to undulate; therefore, we need to consider wave- 
induced modulation of the local direction vectors, s and k. Let the instantaneous position 
vector of the loop be described as a time-varying perturbation about a static equilibrium. 
We will consider two transverse displacements of the axis, one displacement £(s, £) polarized 
in the direction of the principle normal fco (which is confined to the x—z plane) and a 
perpendicular displacement C( s )^) polarized in the binormal y direction (i.e., horizontal and 
parallel to the photospheric plane), 

r(s, t) = r {s) + £(s, t) k {s) + ({s, t) y . (4.22) 

In the equation above, as before, the coordinate s is the pathlength along the unper- 
turbed position of the tube. If we carefully consider the perturbations to the unit vectors, 
we find to first order in the displacements 



s + — k + — y , (4.23) 



= + • (4-25) 
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Since, the equilibrium configuration of the loop lacks curvature in the y direction, the radius 
of curvature does not depend on the ( displacement, at least to linear order. Because of 
this, the two polarizations decouple, and £ and ( satisfy independent equations. These 
two wave equations are obtained by noting that transverse oscillations of the thin fibril 
are largely incompressive; therefore, the longitudinal component of the equation of motion, 
equation (12. ip . can be ignored. 

Differential equations describing the two polarizations of oscillation can be derived by 
inserting equations (I4.23p - (l4.25p into equation ( 12. ip . Fourier transforming in time, and lin- 
earizing in the displacements, 



2 ( d2 1 \ c n 



The influence of the geometry of the equilibrium loop model is felt only through the radius 
of curvature Rq and through the parallel component of gravity g\\ = g ■ s . The terms 
proportional to Cg represent magnetic tension, with the radius of curvature corresponding 
to the contribution from the curvature of the equilibrium field lines. The terms with g\\ are 
buoyancy forces and the inertial terms are those that depend on the temporal frequency u. 
These are quite general equations that describe oscillations for any loop equilibrium — not 
just the isothermal model we consi der here. In fact, excep t for the inclusion of curvature, 



these are the same equations used by I Jain fc Hindmanl ( 120121 ) to study the effects of buoyancy 



and wave speed variations on mode frequencies. 

Let's first examine the tension terms in detail. When these two equations are compared, 
we see that the curvature of the tube only appears in the equation for the normal oscillations. 
This additional tension term from curvature typically has the opposite sign from the tension 
arising from undulations. Thus, the two tension forces that appear in the equation for normal 
oscillations act in opposition and we should expect that curvature of the loop reduces the 
eigenfrequencies. Further, if we consider the limit of small Magnetic Bond Number e, we 
can see from equation ( I3.16P that the radius of curvature is large R = 0(e _1 ). Therefore, 
the curvature of the tube should generate a second-order frequency shift Aw/w = 0(e 2 ). 

If we insert g» = —gz ■ s into equation (I4.27p . and include the Magnetic Bond Number 
dependence directly, we get 
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d 2 C 



ds 2 ~^ X ~° ds~^~ Cn^ 



(4.2* 



A cursory examination of this equation would seem to indicate that for small e, the frequency 
shift from buoyancy forces should be a first-order quantity. However, this is incorrect. We 
have previously demonstrated that the equilibrium height and shape of the loop sensitively 
depend on the Magnetic Bond Number, and in particular as e — > the loop becomes flat 
and horizontal. In fact, by using equations ( 13 .4p . (13. 7p . and ( I3.14p one finds 



z ■ s = sin (ex/ X) . 



(4.29) 



Since sin (ex/A) = 0(e) for small e, it is now clear that buoyancy should cause frequency 
shifts that are second order in the Magnetic Bond Number. This makes sense, as curvature 
also enters at second order and buoyancy and curvature are balanced in the equilibrium 
configuration of the loop. 

The quantities g\\ and R depend on the equilibrium shape of the loop, which depends 
explicitly on the horizontal coordinate x — see equations ( I3.13p - ()3.17p . Therefore, it is quite 
natural to make a change of variable from the arclength variable s to the horizontal x 
coordinate. After inserting the properties of the equilibrium, we find that the resulting 
ODEs lack first-derivative terms, 



d 2 ^ 
dx 2 



+ 



10 



+ 



c 2 cos 2 (ex/A) A 



dx 2 



+ 



cos 2 (ex/A) 



. 



(4.30) 
(4.31) 



The absence of the first derivatives occurs because the tension and buoyancy forces in the 
equilibrium are in exact balance. 

Both of these equations have solutions in the form of Associated Legendre Functions, 



X 



I COS I 



.1/2 



AiPf^ (sin0) + A 2 Q>; /2 (sin i 



C(x) 



cos ( 



,1/2 



A 3 P^ 1/2 (sm6) + A 4 Q II _ 1/2 (sini 



(4.32) 
(4.33) 



with 



0{x) = e|, (4.34) 

/ i , ,2 r 2 \ 1/2 

The dependence on the frequency is buried in the upper index \i of the Associated Legendre 
Functions. Therefore, these solutions are not particularly friendly, and applying boundary 
conditions to quantize the eigenfrequencies would require a nontrivial numerical computa- 
tion. We avoid this difficulty by realizing that the Magnetic Bond Number e is a small 
dimensionless parameter, thus allowing the use of perturbation theory. 

The equations that describe oscillations in both the normal and binormal directions have 
very similar form. In fact, a brief inspection of equations (I4.30P and (14.3ip reveals that the 
eigenf unctions of the two polarizations should be identical up though second order in e. This 
means that the derivation of the solutions for the two polarizations are similar; therefore, we 
will demonstrate the perturbation analysis for only the normal component and simply state 
the result for the binormal component. Assume that the eigenfunction and eigenfrequency 
can be written as perturbation expansions, 



£ n (x) = ~ n (x) + e 2 6Ux) + 0(e 4 ) , (4.36) 

u£ = ^ + 6 2 &4 + 0(e 4 ) . (4.37) 

Insert these expansions into equation (I4.30p and solve in the usual fashion for the zeroth- 
order eigenfunction and eigenfrequency, noting that the cosine appearing in the second term 
can be expanded for small argument, 



E n (x) = Xsm[K n (X-x)\ , (4.38) 

TITCCq 

ll n = n n c = —-r . (4.39) 
ZJv 

The boundary condition of vanishing displacement at the footpoints £(±X) = has been 
applied. 
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If we now collect terms of the next order, we find the following inhomogeneous equation 
for the perturbed eigenfunction, 



dx' x 



+ 



+ 



5ul 



Cl X 2 C 



+ 



X< 



(4.40) 



The particular solution can be found by judicious guessing and by enforcing that the per- 
turbed eigenfunction has the same symmetry about the center of the loop as the unperturbed 
eigenfunction, 



X' 



^n(x) = - — E n (x) 



X 



Co x (5u 2 n 1 



6X 2 ill 2 V C 







2X 2 



KM ■ 



(4.4i; 



The perturbed eigefrequency is fixed by reinforcing that the eigenfunction vanishes at the 
footpoints for all orders of expansion 5£ n (±X) = 0, 



5u„ = - 



Q 2 r 2 
3 2X 2 



(4.42) 



If we now put all the pieces together, to second order in the small parameter e, the 
eigensolution for the normal displacement is given by 



^r,(X 



2 2 
K n C 



l-e 2 



9 / 9 

tX-j __ j iXj I tXj 

JO^ n[x)+ Q V X 2 ' 



0(e 4 



+ 0(e 4 



2<X 2 



(4.43) 
(4.44) 



By the same procedure the eigensolution for the binormal polarization can be obtained, 



Ux) = Ux)+0{e A ) 



2 2 
K n C 



1-e 2 



2n 2 n X 2 



+ 0(e 4 



(4.45) 
(4.46) 



Figure |3] illustrates the resulting eigenfrequencies as a function of e for displacements in both 
the normal and binormal directions. Figure H] shows sample eigenfunctions. 
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5. 



Discussion 



In the previous sections we have modeled a coronal loop as an isolated, thin magnetic 
flux tube. The effects of gravity and curvature have been self-consistently included, leading 
to a static equilibrium determined by the balance of magnetic buoyancy and tension. We 
have computed the eigensolutions for transverse oscillations of both polarizations. In the 
following subsections we will discuss in detail the effects that gravity and curvature have on 
the structure of the eigenfrequencies, assess the differences between the two polarizations of 
wave motion, and explore the implications of relaxing several of our operating assumptions. 

5.1. The Effects of Gravity and Curvature on the Eigenfrequencies 

There are only two free parameters in our model, the footpoint separation 2X and 
the Magnetic Bond Number e. The footpoint separation just provides a physical scale for 
the problem and the structure of the solutions do not depend upon it. The solutions do 
depend intrinsically on the Magnetic Bond Number, because the equilibrium shape of the 
loop is sensitive to the strength of buoyancy and magnetic tension. As the magnitude of 
the parameter e becomes larger (more negative), buoyancy forces increase in strength, which 
causes the loop to rise and a corresponding increase in the curvature of the loop is needed to 
allow magnetic tension to bring the loop back into equilibrium. This results in three ways 
by which the frequencies are modified: the loop lengthens, the restoring force of buoyancy 
increases, and the increased curvature modifies magnetic tension. We can identify the first 
of these effects on the eigenfrequencies by expanding the loop length L, equation (|3.18p . for 
small e, 



If the loop were straight and without gravitational effects, the eigenfrequencies would be 



As expected, an increase in the loop length (increasing e 2 ) decreases the eigenfrequencies 
and the resulting fractional frequency shift is dispersionless and independent of mode order, 




(5.47) 





(5.49) 




3 
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This effect appears directly, in equations (I4.44p and (I4.46P as the first term in the parentheses 
on the right-hand sides. 

Buoyancy appears in identical form in the equations describing each polarization of 
transverse oscillation; both equation ( 14 . 2 6 [) and (I4.27P have the same first- derivative term. 
This common form, leads to identical perturbations to the eigenfrequencies due to the in- 
clusion of buoyancy. From equations (14.441) and (14.461) we can directly see that buoyancy 
increases the eigenfrequencies by the fractional amount, 



^ - e 2 ^- (5 50) 

This might be a somewhat surprising result. A simple dimensional analysis of the equation 
of motion would indicate that buoyancy causes frequency shifts t hat are first orde r in th e 



Magnetic Bond Number e, and indeed this is what was found in iJain fc Hindmanl (120121 ). 
However, the buoyancy force in the wave equation is proportional to e times the geometrical 
factor z- s . This geometrical factor also depends on e because the shape of the coronal loop 
is sensitive to e. Low Magnetic Bond Numbers correspond to loops that are low-lying and 
nearly horizontal, and because of this the geometrical factor is first order in the Magnetic 
Bond Number, and the buoyancy force in total is second order. 

The last effect, curvature of the loop, modifies the magnetic tension than can act as 
a restoring force. It only appears in the equation for the normal polarization, and its sign 
indicates that it counteracts the local curvature induced by the wave motions (i.e., the second- 
derivative of the displacement). Therefore, we would expect curvature of the background 
loop model to decrease the eigenfrequencies. This does indeed occur, as we can see directly 
from the difference between the frequencies of the normal and binormal polarizations. The 
curvature of the loop results in a negative frequency shift that is twice as large in magnitude 
as the effect of buoyancy, 



# - -*4? ■ (5-51) 

The frequency shifts induced by both buoyancy and curvature enter at second order in 
the Magnetic Bond Number. That both enter at the same order is a simple consequence 
of the balance between buoyancy and tension in the equilibrium. That they both enter at 
second order could easily be predicted. From dimensional arguments alone, the shift in the 
square frequency due to curvature must scale with the radius of curvature as Acu 2 ~ c^Rq 2 . 
Since the radius of curvature is inversely proportional to e, the frequency shift must be second 
order in e. 
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A consistent result for the curvature was found by Ivan Doorsselaere et al.l (120041 ) where 
they solved for fast waves in a toroidal coordinate system in the absence of gravity. The 
equilibrium field was assumed to be purely toroidal and a coronal loop modeled by a partial 
section of the torus. With an extensive mathematical analysis they concluded that the 
curvature of their torus induces a frequency shift that vanishes at first order in the ratio of 
a/Ro, where a is the minor radius of the torus (and Ro corresponds the the major radius). 
Therefore, the first non-zero perturbation must occur at second or higher order. Here we 
have shown that for our thin tube model the first non-zero perturbation to the frequency 
does indeed occur at second order in the radius of curvature, as a simple dimensional analysis 
would indicate. 



5.2. The Two Polarizations of Oscillation 



Many of the observed transverse oscillations appear to be horizontal (e.g.. lAschwanden et al 



20021 ). and therefore in the binormal direction. Unfortunately, our assumption of an isolated 
magnetic loop is particularly suspect for such oscillations. Most of what we call coronal loops 
are probably just a bundle of field lines that are part of a much larger magnetic structure like 
a magnetic arcade. We identify this particular bundle as a "coronal loop" solely because it 
is selectively illuminated by a local heating process. If the loop is part of a magnetic arcade, 
oscillations in the direction of the principle normal (sometimes called "vertical" oscillations) 
may be unaffected by the larger arcade structure, because the undulations do not impinge on 
neighboring arcade field lines. However, binormal (horizontal) oscillations must inherently 
impact nearby field lines and the anisotropic nature of the arcade structure must modify 
wave propagation. In particular, waves may no longer be confined to the loop. Spanwise 
propagation along the arcade is possible, and since the other loops in the arcade may not be 
illuminated, this propagation may go undetected. Such waves may still be trapped, but the 
problem is inherently two-dimensional and our one-dimensional analysis incorrect. 

Despite these difficulties, there are distinct differences between the two polarizations 
of waves that may be observable. Unfortunately, the eigenfunctions are of little direct use 
because the two polarizations have identical eigenfunctions (to second order in the square 
of the Magnetic Bond Number). Therefore, observations of loop motion must be able to 
directly account for projection and occultation effects if they are to determine polarization. 
The frequencies of oscillation are not degenerate, however, and offer promise to winnow 
between polarization modes. 



Many observational and theoretical st udies have examined the diagnostic value of the 
ratio of eigenfrequencies (for a review see lAndries et al.l 120091 ) . in particular the ratio of 
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the fundamental and first overtone. For a straight magnetic fibril in the absence of gravity 
and stratification, the ratio of the first overtone to the fundamental U2/001 is 2. Observa- 
tional differences fro m 2 have been interpreted as evidence for wave speed stratification (e.g., 



McEwan et al.ll2008[ ). and potential wave speed models have been rejected based on whether 



the model generates a ratio tha is less than or greater than 2. As we shall soon see, curvature 
and buoyancy forces can also cause the ratio to shift from the canonical value of 2, without 
any spatial variation in the wave speed. 

For our model, the frequency ratio for any given pair of neighboring modes is given by 



(5.52) 

where the upper sign (+) corresponds to oscillations in the normal direction and the lower 
sign (— ) to those in the binormal direction. For the ratio of the first overtone to the funda- 
mental, this becomes 



n 



1 ± 



n 



e 2 2n + l 
7T 2 n 2 (n + iy 



+ 0(e 4 



- = 2±|^ + 0(e 4 ). (5.53) 

Since the effects of the variation in the loop length are dispersionless, they have dropped 
out in these ratio expressions, leaving just the effects of buoyancy and curvature. Direct 
buoyant augmentation of the restoring force causes the ratio to decrease, while curvature 
of the loop causes the ratio to increase. Since the normal polarization senses the effects of 
loop curvature, whereas the binormal polarization does not, the ratio of the first overtone 
to the fundamental is greater than 2 for normal oscillations and less than 2 for binormal 
displacements. The diagnostic value of the frequency ratio is clearly predicated on the 
knowledge of which polarization one is observing. 



5.3. The Effects of an External Magnetic Field 

While we have assumed that the coronal loop is embedded in a field-free atmosphere, 
we can speculate what the effects of an external magnetic field might be. Let's first consider 
the equilibrium. The inclusion of an external magnetic field only provides minor modifica- 
tions to the steady state form of the equation of motion ( 12. ip . The tangential component, 
equation (13. 8p . remains unchanged as it depends only on internal fluid properties. The same 
can be argued for the magnetic tension as long as the thin-tube approximation holds. The 
buoyancy force is modified, however, because the density difference between the internal and 
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external fluid is a function of the external magnetic field through equilibration of the total 
pressure across the tube. However, if the external field is force-free, the density scale height 
of the external fluid is unaffected and density changes only manifest as a different value of 
the Magnetic Bond Number e — see equation (13 .lip . Therefore, the equation that describes 
the equilibrium shape of the loop, equation (I3.10p . remains valid in the presence of external 
field, and of course the solutions of that equation also remain valid. 

The existence of an external magnetic field is not as benign for the wave equations 
satisfied by the kink waves. The presence of the surrounding magnetic field can be felt by 
the waves through backreaction forces caused by the tube's motion. The well-known effect 
of enhanced inertia, which accounts for the fact that a moving tube must move external 
fluid out of its way, has already been included in our equation of motion (I2.ip . However, 
when the external fluid is magnetized, there is an additional effect because bending of the 
tube requires bending of the external field as well. Therefore, the external field B e provides 
additional stiffness that increases the kink wave speed, 



This well-known expression for the kink wave speed (e.g.. lEdwin fc Robertslll983l ) displays 



both backreaction forces explicitly. In the denominator the sum of the internal and external 
densities is the effect of the enhanced inertia, while in the numerator the appearance of B 2 
is the result of the increases stiffness. These forces are derived and identified explicitly in 
Appendix [A] for a thin tube in the absence of gravity. 

When including the effects of an external magnetic field, one is tempted to simply 
replace the kink speed Co appearing in the wave equations (I4.26[) and (14.271) with the proper 
wave speed ck, 



d 2 j gjDydt 
ds 2 c 2 K ds 



+ 



uj* 1 
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ds 2 c 2 js ds ' c 2 



(5.55) 
(5.56) 



From the motivating calculation in Appendix |A] it is evident that such a substitution is 
justified in the absence of curvature of the loop. When curvature is present, the validity of 
the substitution seems physically appropriate, but remains mathematically unproven. While 
these two wave equations possess the same form as those solved previously, equations (I4.26P 
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and (14.271) . there is one subtle difference. The ratio appearing in front of the first derivative 
in each equation is no longer related to the parameter e in the same fashion. Previously, the 
denominator depended on cq instead of Cr-. Since, the magnetic tension in the equilibrium 
still depends on c (and not cr-), after inserting g\\ appropriate for the equilibirum, we find 



,2 



-K ±L 

d 2 ( t iv\dC u 



^ + a - sm{ex ,X)f s + l- r + M ) ( = 0, (5.57) 



^ + = 0, (5.58) 

with a being a factor that depends on the plasma parameter within the tube, /3, and outside 
the tube, (3 e , 



a = WTbI = JTWTK ■ (5 ' 59) 

Since the corona is magnetically dominated (i.e., (3 « 1 and (3 e « 1), a is generally quite 
close to a value of one-half. The form of our previous wave equations can be recovered only 
for a = 1. While this seems like a minor difference, the change of variable that we performed 
to remove the first derivatives in equations (14.301) and (I4.3ip is predicated on a = 1. So, 
while the basic form of the wave equations remains unchanged when we include an external 
magnetic field, the detailed perturbation analysis performed in §3]is n o longer valid. It would 
be relatively trivial to repeat the perturbation analysis for this new case, but we choose to 
leave such an effort to a future study. 



6. Conclusions 



We have developed a model of a coronal loop that self-consistently includes the effects 
of gravity and curvature. This has been accomplished for a loop that is treated as an 
isolated magnetic fibril embedded in an isothermal field-free corona. For such a model, the 
shape adopted by the loop when in a state of static equilibrium is completely determined 
by a balance between the forces of buoyancy and magnetic tension. The height of the loop 
depends on the Magnetic Bond Number e which is usually small for coronal loops. For 
these small values, the loops are low-lying with weak curvature. In the other extreme, no 
stable equilibrium exists for |e| > n/2. This stability criterion is actually a restriction on the 
footpoint separation compared to the coronal density scale height. Loops with a footpoint 
separation greater than 2irH are unstable. 
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Since the static loop is confined within a vertical plane, and only possesses curvature 
within this plane, the two polarizations of transverse oscillations (those in the direction of the 
principle normal and those in the direction of the binormal) decouple. The wave equations 
that describe these oscillations both have restoring forces from buoyancy and from the mag- 
netic tension induced by undulations. The equation that describes the normal polarization 
also possesses a tension term that arises from the curvature of the loop itself. The effects of 
the restoring force of buoyancy and of background curvature both enter the eigenfrequen- 
cies at second order in the Magnetic Bond Number. Buoyancy increases the frequency and 
curvature reduces it. The ratio of the frequencies of two eigenmodes of neighboring order 
u n+ i/u n differ from the canonical value of (n + l)/n that would be achieved for a straight 
loop in the absence of gravity with constant kink speed. Oscillations confined to the plane 
of the loop, or those polarized in the normal direction, have a frequency ratio that exceeds 
this canonical value, u n+ i/u n > (n + l)/n, whereas horizontal oscillations, or the binormal 
polarization, have a ratio that is less than the canonical value, u n+ i/u n < (n + l)/n. These 
conditions on the frequency ratios are achieved for a loop with a constant kink wave speed, 
independent of the pressure scale height. 
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A. Backreaction Forces on a Magnetized Tube in a Magnetized Environment 



In order to see how the backreaction forces might appear for a magnetic fibril moving 
through a magnetized external fluid, we perform a simple wave emission calculation for a 
field configuratio n amenable to analytical t reatment and direct identification of the relevant 
forces. Following Edwin &: Roberts! (119831 ). consider a magnetized cylinder embedded in a 
magnetized atmosphere in the absence of gravity. The magnetic field within the cylinder is 
constant and points in the z-direction, B = Bz. The magnetic field in the region exterior 
to the cylinder is also constant and points in the same direction B e = B e z. While the two 
fields are aligned, they may potentially differ in strength. The mass densities inside and 
outside the cylinder are p and p e , respectively. 

Magnetosonic waves outside the cylinder have a simple analytic form when expressed in 
polar coordinates with the axis of the coordinate system coaligned with the axis of the cylin- 
der. A general for m for the radial fluid di splacement £ r and the total pressure perturbation 
611 are as follows (jEdwin fc Roberts! Il983l ): 
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£ r = k- 1 ?™* Q' m (Ar) e iKZ e~ iuJt , (Al) 

SU = p e (Xny 1 (tu 2 - k 2 V 2 ) e im * Q m (Ar) e iKZ e~ luJt , (A2) 
with the dispersion relationship 

_ (a, 2 - „ 2 c 2 ) (u 2 - k 2 V 2 ) 

(Cl + V 2 ) (W 2 - K 2 U*) ■ 1 1 

In these equations, c e , V e and U e correspond to the external values of the sound, Alfven, 
and slow wave speeds, respectively. The axial and radial wavenumbers are denoted k and A, 
and the temporal frequency is u. The function Q m is a mth-order Bessel function or Hankel 
function. Primes denote differentiation with respect to the argument of the Bessel function. 

If an axially propagating wave causes the tube to oscillate back and forth in the x 
direction with an amplitude rj, 

£tube = V e iKZ e~**& , (A4) 
the radial displacement at the lateral surface of the tube (r = a) is given by 



er,tube = r/cos0e^e-^. (A5) 

The radial displacement of the external fluid must match the normal motion of the tube's 
outer surface. Therefore, the wavefield in the external medium must have the same cos0 
azimuthal dependence as appears in equation (1A5j) and will be composed of only the dipole 
components (m = ±1). Further if A 2 > 0, causality requires that the wavefield only possesses 
outward propagating waves, leading to the following form: 



£ r = Ak~ x cos 4> H[ (Ar) e iKZ e~ iut , (A6) 

5U = An- 1 Pe A" 1 (u 2 - k 2 V 2 ) cos0 H x (Ar) e iKZ e~ iuJt , (A7) 

with the function Hi corresponding to an outward propagating Hankel function of the first 
kind, Hi(Xr) = ^ (1) (Ar), and the wave amplitude A fixed by the displacement matching 
condition at the tube's surface, 
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The backreaction forces acting on the tube due to the tube's motion through the external 
fluid can be derived by integrating the forces acting on the tube over the tube's cross-section. 
In the absence of gravity the only forces acting on the tube are magnetic tension and total 
pressure. If the tube is thin, we may ignore internal variations of the magnetic field strength 
B and the cross- sectional mean of the tension force acting on the tube is given by 



tension 



D a Stube _ 

47r dz 2 

where V is the Alfven speed within the tube. 



- P V 2 K 2 r] e iKZ e 



-iojt 
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(A9) 



The cross-sectional mean of the total pressure force can be expressed as the integral 
of the total pressure perturbation around the circumference of the tube by using Stokes 
theorem to convert the area integral into a contour integral, 



pressure 
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The x-component of the total pressure force is the only nonvanishing component, and the 
mean force in this direction is therefore given by the sum of the pressure and tension forces, 



= -- 



2- 



cos (j) oTI(a, 0, z) d(f) — pV k rj e 



2,„2„ „ikz „—uct 



(AH) 



The integral appearing in this equation is trivial because the total pressure field on the 
surface of the tube only has dipole dependence. Using equations f]A7|) and (IA80 this reduces 
to 



F„ = — 



P kV r] + p e {uj -kV*) (Aa) jph^i V 



(A12) 



Since the tube is thin we can perform small argument expansions on the Hankel functions 
to obtain a simple form, 
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K 4^ 



p e uj 
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(A13) 
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The equation of motion can therefore be written 



J^2 _|_ ^2 

- P^tube = 4^ e rf Z 2 ^tubc + PeW 2 ^ tube . (A14) 

We can immediately see that there are two backreaction forces. Both arise from the 
total pressure in the external fluid. The term proportional to p e w 2 £ tu be i s the familiar effect 
of enhanced inertia and arises from the u 2 factor appearing in equation (1A7|) . In the absence 



of external field, this would be the only term. However, when external field is present, 
the total pressure in the external fluid is reduced by the fractional amount (w 2 — k 2 V 2 ) /u 2 
compared to pressure in the unmagnetized fluid. The term proportional to the square of the 
Alfven speed leads to an additional backreaction term — K 2 (£> 2 /47r)£ tubo that represents the 
force needed to bend the external field lines when the tube itself bends. The external field 
provides additional stiffness that increases the wave speed. 
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Fig. 1. — The geometry of the coronal loop, with the loop shown in black. The photosphere 
corresponds to the x—y plane, while the loop is confined to the x—z plane. The local Frenet 
coordinates are indicated in red. The local tangent vector is sq and the principle normal fco lies in 
the direction of curvature. The binormal is everywhere constant and pointed in the y direction. 
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Fig. 2. — Equilibrium (a) height and (b) radius of curvature of coronal loop models with differing 
ratios of buoyancy to magnetic tension, i.e., differing values of the Magnetic Bond Number e. The 
value of e associated with each color is indicated in the right panel. In both panels the dashed line 
corresponds to a semi-circle with constant radius of curvature Rq = X. Loops with weak buoyancy 
(|e| << 1) are flat with large radius of curvature everywhere. Loops with substantial buoyancy 
(e ~ — 7r/2) are tall with large radius of curvature in the legs and small radius of curvature at its 
apex. Loops with positive e or with e < —it/ 2 are unstable. 
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Fig. 3. — The eigenfrequencies of coronal loop waves as a function of the Magnetic Bond Num- 
ber e. (a) Fractional perturbation in the square eigenfrequency for both polarizations of motion. 
Transverse oscillations in the normal direction are shown with the solid curves, while those in the 
binormal direction are illustrated with dashed curves. The different colors correspond to different 
mode orders: the fundamental (n = 1) shown in black and the overtones in red (n = 2), green 
(n = 3), blue (n = 4), and aqua (n = 5). The dotted curve corresponds to the fractional frequency 
shift (independent of mode order) that arises solely from changes in the length of the loop, (b) 
Frequency ratios of nearby mode orders, with the same line styles and colors as above (i.e., black 
corresponds to u^/oJi). The dotted curves indicate the value appropriate for a straight loop in the 
absence of gravity (n + l)/n. 
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Fig. 4. — Eigenfunctions for the fundamental (black) and first two overtones (red and blue) as 
a function of arclength along the loop. The solid curves illustrate the waveform for a Magnetic 
Bond Number of e = — 7r/4, whereas the dashed curves correspond to those for the unperturbed 
eigenfunction with e = 0. Both displacement polarizations, £ and possess the same eigenfunctions. 




